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ABSTRACT 

We show that partition functions of various matrix models can be obtained by acting on elementary functions with 
exponents of W operators. A number of illustrations is given, including the Gaussian Hermitian matrix model, Hermitian 
model in external field and the Hurwitz-Kontsevitch model, for which we suggest an elegant matrix-model representation. In 
all these examples, the relevant W operators belong to the iy( 3 ) algebra. 

Contents 

1 Introduction 1 

2 Hermitian Matrix Model in external field and operator W2 2 

3 Gaussian Hermitian Model and operator W 2 5 

3.1 The main relation, eq.(4) 5 

3.2 Eq. (4) from Faddeev-Popov trick 7 

3.3 Eq.(4) from Virasoro constraints 8 

3.4 Non-Gaussian models and operators W- p 9 

4 Hurwitz-Kontsevich Model and operator Wq 9 

4.1 Approximate matrix integral 10 

4.2 Corrections 11 

4.3 Exact matrix integral 13 

4.4 Character expansion 16 

4.5 Appendix: GL(n) characters as eigenfunctions of Wq 18 

5 Conclusion 22 



1 Introduction 

These days we witness a renaissance of matrix model theory, both of its applications and of "theoretical 
theory" . Matrix models are finally recognized as a source of new, badly needed, special functions and as a 
simplified, still representative, model of entire string/M-theory, with sophisticated generalized geometries 
of Calabi-Yau type behind the vacua structure substituted by a far better studied geometry of the Ricmann 
surfaces. For summaries of the previous stages of matrix model theory see, for example, [l]-[5] and references 
therein. For recent papers with a number of advanced new developments see [6] -[17]. 

One of important issues about string theory partition functions is their generation by canonical proce- 
dures from simple canonical objects. There is a whole variety of such reductionistic properties. In the case 
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of matrix models one can think of reducing them to a simpler (more fundamental) theory - like that of 
free fields on Ricmann surfaces and thus to 2d conformal and finally to group theory, see [2, 5, 10, 11] for 
different stages of this project. One can embed matrix models into group theory in a somewhat different 
way, by exploring and exploiting integrability properties of partition functions [3, 5] (partition function of 
a quantum theory is always a kind of a r-function [18], and in the case of the eigenvalue matrix models 
these are usually r-functions from the well studied KP-Toda family, associated with the U(l) Lie algebra). 
One can instead express generic matrix models through a few simple ones, like the Kontsevich model [3], 
which has alternative origins in combinatorics and geometry of moduli spaces, see [13]- [15]. This is a part 
of reductionist program within the matrix model field itself, especially important if one uses it to model 
the pertinent features of string theory [13]. 

In this paper we address another reduction of the same type: from sophisticated to simple r-functions, 
but a much simpler one as compared to meron/instanton decompositions of [13]. Namely, as previously 
observed in the case of the Hurwitz-Kontsevich model [14], partition functions can be generated from some 
trivial r-functions - like e* 1 - by the action of non-trivial generators 

Partition Function = e w ( Elementary Function ) 

from integrability-preserving GL(oo) group, which converts one family of Virasoro-like constraints into 
another (it is actually enough to look at the string equations). See the basics of the underlying theory of 
equivalent hierarchies in [19]. Operators W are naturally classified by their spin: when constructed from 
free fields, the spin-fc operators have the form 

m— — oo 

i.e. are made from the fc-th powers of the U(l) currents on Ricmann surfaces. In the simplest examples, 
which we consider in this paper, the relevant generators are just the next-complicated after the spin-2 
Virasoro ones: the operators [20] (we call them simply W in what follows). When expressed through 
the n x n matrix- valued background field ip (the Miwa variable), operators are differential operators 

of order (k — 1), so for k = 3 they resemble Laplace operators: 



W 2 (3) = tr 



d V o 2 



and 



dip J dtydipl 



dip* dip : 



are the simplest illustrations of this property, which plays an important role below. The goal of this paper 
is just to describe a few examples, leaving intriguing applications to separate publications. Some relations 
to GKM theory [3] are immediately obvious, but even they will be discussed elsewhere - in order to clearly 
separate explicit formulas of the present paper from broader hypotheses and speculations. 

2 Hermitian Matrix Model in external field and operator W 2 

To begin with, we consider external-field correlators in the Hermitian matrix model 



Ck u ...,k m - 

NxN 



J dep 4,2/2 Tr(cp + tp) kl ...Tr(<p + tP) k 
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where the external field ip is a constant Hermitian N x N matrix and d<fr — j d4> % j is the flat measure. 
Obviously, correllators C kl ,..., km are invariant under conjugation ip UipU^ 1 , i.e. they are functions of 
invariant variables T k = Trtp k . For example 

C x =T a 

C 2 =T 2 + T 2 

C 3 = T 3 + 3ToT! 

Ci,i = ?i 2 + T 

C 2 , 2 = T| + 4T 2 + 2T 2 T 2 + 2T 2 + T 4 

Ci,i,2 = T 2 T 2 + T T 2 + 4T 2 + T 2 T 2 + 2T + T 3 

Note, that correlators depend on N only through T = Tr tp° — N. With the help of the shift operator 

exp(Tr<^!^/(V>) =/(<£ + V) V/, 

where (tp T Y- — ipl is the transposed matrix, the integral over can be made Gaussian: 



C kl ,...,k m = j 'd<f> e - Tr ^ 2 /2 Tr(^ + V) fel ...Tr(^ + V) fem 

= y d0 exp ^-Tr 2 /2 + Tr Tr . . . Tr ^ 



° XP (V (^) 2 ) Tr ^ fel ■ • - Tr ^ m = ex P (V (^)l Tr ^ • ■ - Tr ^" 



In this way we find an explicit formula for all correlators: 



C kl ,...,k m = e w '' 2 T kl ...T km (1) 

where operator 

d \ 2 d 2 



W 2 =Tr 



,3 



can be called a matrix Laplace operator. It converts invariant (under conjugation of ip) functions into in- 
variant functions, and therefore, can be reduced to the space of such functions, where it acts as a differential 
operator of second order in invariant variables T k . Indeed, by application of the chain rule 

d p(T) V dT - dF ( y ) 
tyi \ > Z^i dTa 
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and similarly 



F (T)= V dTa dT " d2F{T) V ^ dFiT) 
a j^ dr jdi> {dT a dT b ^itydti 9T a 



Taking derivatives of traces, it is easy to check that 



dT a 8T b 



T = abT a+b - 2 



and 



Therefore 



d 2 T a 



{k + l + 2)T k Ti 



3 Y i k+l=a-2 



(f) \ 2 00 / Ci 
7TT = E [(a + b + 2)T a T b — + abT a+b - 2 



d 2 



dT a 8T b 



(2) 



Once again, the last identity is true, when the operator acts on invariant functions, i.e. on functions of 
time- variables T k . As usual, the partition function is introduced as generating function for all correlators. 
It depends on two sets of time- variables, t k and T k : 2 

Z W = E h. E C kl ,... tkm (T) t kl ...t km = d<j> e-^+S^W^ 

m=0 fei,...,fe m =0 

It follows from (1), that 



Z(t\T) 



,W 2 /2 p J2 k>0 t k T k 



(3) 



As one can see, partition function of the Hermitian matrix model in external field is generated from the 
trivial function e^ tkTk by the action of generator W 2 /2. Formula (3) is quite interesting: an explicit 
representation for the partition function, which does not involve matrix integrals. Also, eqs. (1) and (2) 
are very convenient to calculate particular correlators "by bare hands". For example, 



G 



2,2 



l + ^W 2 + -{W 2 ) 2 )Tl 



since all powers of W 2l higher than 2, annihilate Tf . One finds 

W 2 T 2 2 = 8T 2 + 4T 2 T 2 , 



(W 2 ) 2 Ti = 16T 2 + 8T 4 



and 



C 2 , 2 = T| + 4T 2 + 2T 2 T 2 + 2T 2 + T 4 



2 Through this section, to simplify our formulas and to make closer contact with [14] in s.4 below, we omit the factor 1/fc in 
Miwa transform T k = Tri/> fe , what, actually, spoils the natural symmetry between f fc and Tfc. If 1/fc is restored, T k = Tr?/) fc /fc, 
then the exponent in (3) acquires its usual form exp (fct fe r fc ). In this case, however, one gets an additional factor exp (to To). 
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3 Gaussian Hermitian Model and operator W-2 
3.1 The main relation, eq.(4) 

Our next example is the simplest matrix model - Hermitian matrix model in Gaussian potential: 

Z G = J d0 e - T ^ 2 /2+2: fc >o* fc Tr^ 
NxN 

where the integral is taken over all N x N Hermitian matrices with flat measure and depends on the set of 
time- variables tk- Since Tr0° = N, the dependence on to is given by simple multiplicative factor e^* : 

Z G = e Nt ° J dtfe- 1 ** 2 / 2 +£*><> 

WxJV 

We are going to demonstrate, that partition function of the Gaussian model is generated by operator W-2 



Z G = e w -^ z e 



W-2 /2 Nt 



(4) 



where 

W- 2 = V abt a t h — + (a + b + 2)t a+b+2 — — (5) 

T^„ V Ota+b-2 dtadtbj 

a, 6=0 

Note, that W2 in (3) is acting on T, while W-2 in (4) - on ^-variables. To prove (4), we make a Miwa 
transform - introduce an n x n Hermitian matrix -0, such that 

tk = T tr V'~ fe , k > 
k 

It is important that the size n of matrix ip is absolutely independent of the initial size N of matrix <p, 
because interaction terms in the action do not involve addition or multiplication of these matrices (this 
was not the case in s.2). This is the usual feature of Kontsevich-like matrix models [3], what emphasizes 
relation between the subject of this paper and GKM theory. Consequently, we distinguish operations tr 
and Tr , which denote traces of matrices of sizes n and N, respectively. 
After the transform, the left hand side of (4) takes form 



J d<i> e -^ 2 /2 exp (j- I tr 0- fe Tr^ j 

w AT \k—l / 



Z G = e" t0 

NxN 

Note, that to is not affected by this transform and remains a free parameter. Using the identity 

(°° 1 \ 
-trV' _fe Tr0 fe J 

the interaction terms are written in determinantal form: 

Z n = e Nto [ dSe~ T ^ 2 ' 2 - = e Nta [ dS e~ T ^ 2 / 2 ^ det ^ (6) 

J det (/ ® I - V'" 1 ® <j>) ' J det U <8> / - / ® <£) 

NxN NxN 
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where (<j) ® i>) % jl = is the tensor product of two matrices, regarded as a nN x nA^ matrix. Note, that 

i, j and n, zv are different sorts of indices: i and j take values 1, . . . , N, while /i and v take values 1, . . . , n. 
Eq. (6) is the expression that we need for the left hand side of (4). 

Now let us calculate the right hand side. Operator W-2 can be expressed as a differential operator of 
the second order in terms of tp and to- However, we can always substitute the derivatives with respect to 
to by factors of N. A straightforward application of the chain rule, similarly to the previous sections, gives 




Its exponential is easy to find: 




Using the identity 

oxp (h T {ip)) = I # cxp H tr ^ 2 +tr (^)) 

* ' nxn 

and the properties of the shift operator, we obtain: 

-p(l*-,)^- - /#e X p(-I.^ + ,,(4))^ 

nxn v ' 

/ \ N 

J \det(iP + 4>) 

nxn x ' 

This is what we get for the right hand side of (4). At first sight, this seems different from (6): 

( \ N / \ N 

[ d4> e^ 2 / 2 , vs. / # e-^ 2 / 2 ( d£tV_\ 

J det V ® / - I ® 4>) J \ det N> + 4> ) 

NxN nxn x ' 

even the integration goes over matrices of different size. However, in fact these two integrals are equal and 
we prove this fact in the following subsection. 
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3.2 Eq.(4) from Faddeev-Popov trick 

We now prove (4) by proving the identity between these integrals: 



/ 

NxN 



det (tp ® J - I <8> </>) 7 



JV 



d0 e 



-tr0 2 /2 



det -0 



AT 



det (?/> + 4>) 



or in a more symmetric form 



/ 



d(j) e 



-Tr0 2 /2_ 



NxN 



det (V> ® / - / i 



/ 



d0 e" 



-tr /2_ 



det <8> 7 + <g> /) 



(7) 



We will use Faddeev-Popov's trick, i.e. representation of det 1 as a Gaussian integral over auxiliary fields: 



/ 



d(j) e 



-Tr0 2 /2_ 



NxN 



det (ip ® / - I ' 



dbdc exp ( - ^Tr^ 2 + - <J£$)c"') 



NxN 



-tr<£ 2 /2_ 



det (J 



+ = / W rfMc cxp ( - \ il4>2 + b ^ + 



Here b^i and c v i = 6*^- are bosonic (since determinant stands in denominator) Faddeev-Popov fields. 
Integrals over <j> are Gaussian. After <p is integrated out, we obtain 



/ 



d<j) e 



-Tr0 2 /2 



NxN 



det I - I « <rt - ./ ^ CXP G + 



/ ^ e " W2/2 det(/J + /^) - / ^ CXP (I h ^ j(Mi + 



i.e. the integrals become the same. Thus, our identity 



/ 



d(f> e 



-Tr0 2 /2_ 



NxN 



det (ip 1 



r- = [ d<j> e- tr<i 



2 /2. 



det ® 7 + ^ ® /) 



is true. Therefore, representation (4) is valid. 

At that point it is worth mentioning, that there are other identities between Gaussian integrals, similar 
to (7). For example, there is a direct analogue of (7), with determinants standing in the numerator: 



J d</>e- Tr02/2 det(V®/-/®^) = J d0e +tr * 2/2 det(V>®/ + 0®/) - J dcf> e +tr * 2/2 det (</> + tp) 



N 



NxN 



This identity is well-known as equivalence [21] between the Gaussian model and logarithmic Kontsevich 
model and is usually proved by orthogonal polynomial techniques [5]. Faddeev-Popov's trick seems to be a 
more economic and elegant way to prove such identities. When determinant stands in the numerator, one 
only needs to consider grassmanian Faddeev-Popov ghosts b^i and c„j, as it is usually done in Yang-Mills 
theory. Remarkably, identity breaks down, if determinants are raised to any other power, different from ±1. 
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3.3 Eq.(4) from Virasoro constraints 

There are different other ways to derive (4). For example, Gaussian partition function satisfies a consistent 
system of linear differential equations called Virasoro constraints 

^-Z G = L b ^Z G = ljTat a --^_ + V J-AIzg, b>l; ^-Z G = NZ G 

and (4) is their direct corollary. Indeed, summing by b from 1 to infinity with weight bt b , we obtain 

DZ G = W- 2 Z G (8) 
Two operators appear in this equality: the degree (dilatation) operator 

D = zZ^aT = L o- N2 (9) 



. 'dt q 



and our familiar W-2-operator 



W-2= V [abtatb— + {a + b + 2)t a+b+2 — — 

Vn V ot a+b - 2 at a dt b 

a, 6=0 

with commutation relation 

DW-2 - W- 2 D = 2W-2 (10) 
Notice, that Z G is graded by the total t-degree: 

oo oo 
s=0 m=0 fei + ...+fe m =s 

Operator D preserves this grading: 

DZ^ = sZ { ^ 

As follows from (8), operator W-2 respects the grading in the following sense: 

fV_ 2 4 s) - ( s + 2)4 s+2) 

and this implies that graded components of Z G are generated, one by one, from the lowest component: 



4 2) = -Ji--izg } 



z G = 4°) + i^) + ^(^ 2 ) 2 4°) + ^(^ 2 ) 3 z^ + ... = e --/ 2 4 0) 

where 4^ obviously equals e^' , so that we derived (4) once again. 



3.4 Non-Gaussian models and operators W_ p 

It is tempting to generalize the above Virasoro derivation of (4) to non-Gaussian partition functions 



Z.vg = f d<P e- lrcp - /p +2. fc >o**«o ; p > 2 

JVxAT 

since Virasoro constraints for these models are not very complicated: 

d ( d d d i 

7^;vg= Vai a - + V ——\Z NG , b>(p-l) 

However, in the non-Gaussian case Virasoro constraints are labeled by b > (p — 1), not by b > 1. Immediate 
consequence of this is that Znq is no longer fixed unambiguously by Virasoro constraints alone: some 
additional requirements should be imposed (see [6] for the best studied Dijkgraaf-Vafa example). The 
technical procedure from s.3.3 is also inapplicable, because it is impossible to sum by b from 1 to infinity. 
Instead, one can sum from (p — 1) to infinity and obtain 

D+Z G = W- P Z G (11) 

where 

b=p— 1 6=0 

and 

W- P = V [abtatb— h (a + b + p)t a+b+p — — 

ab^o V dt a+b - p y dt a dt b 

Operator D is the degree operator (9) in all t-variables and satisfies 

DW- P - W-pD = pW-p (12) 

while D + does not satisfy (12) and has the meaning of degree operator in variables tj with i > (p— 2) - only 
a part of all variables. Such a partial grading is not very useful, since W- p does not respect this grading. 
For these reasons, we can not deal with (11) as we did with (8). Some additional ideas are required to 
obtain a e w representation for non-Gaussian, in particular, the Dijkgraaf-Vafa partition functions. 

4 Hurwitz-Kontsevich Model and operator Wq 

In the previous sections, we converted a matrix integral into exponent of a VF-operator, acting on a simple 
function. In this section, an inverse problem is considered: namely, conversion of the e w formula into 
a matrix integral. We discuss this topic with the example of Hurwitz-Kontsevich partition function [14], 
since historically the e w -representation for this function was found before the matrix integral. In result, 
we obtain an interesting matrix model representation for Z^x- 

According to [14], Zhk depends on the time- variables pk and additional deformation parameter t via 



(p) = e tWo ' 2 e^ 



(13) 



where 

OO / 

W n = "V ( („.4-h\r>-n h , . . 

OPa+b OpaOp, 



w =J2(( a + h ^T^r- + »^+& ) ( 14 ) 



a, 6=1 
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Existence of such a formula suggests to look for a matrix integral, responsible for the appearance of Wq. 
We find this integral in several steps: first, we find an approximate matrix integral, then calculate a few 
corrections and finally conjecture an exact matrix integral. In the spirit of the previous example, we 
introduce a Miwa variable - an n x n matrix tp, such that 



Pk = tr^ fc 



Note, that conventionally Miwa transform is defined as tk = tr tp k /k, but we use the rescaled times pk — ktk 
as in [14]. Operator Wo can be expressed as a differential operator of second order in terms of tp. Using 
the chain rule, just like in section 2, we obtain: 

d 2 = y dp a dp b d 2 F(p) | ~ d 2 p a OF(p) 

d^ijdipki dp a dp b f-f dtp i: jdtp k i dp a 

and operator (14) is reproduced if we contract this relation with tp k jtpu: 



Wo = ipkjipil- 



d 2 



oipijdipki 

Let us remind, as usual, that identity 



ntr I tp 



d 



dip 1 



(15) 



d 2 



dipijdtpki 



d d 2 
^ ({a + b)p a pb— Vabp a+b 



a,b=l 



dp a 



+b 



dpadpi 



is true, when the operator acts on invariant functions, i.e. on functions of time- variables pk- 
4.1 Approximate matrix integral 

Having a differential operator in terms of tp, one can rewrite its exponent as a matrix integral over auxiliary 
matrix (p. Unfortunately, it is not possible to literally apply the method of the previous section - the identity 
one would use for this purpose is not quite true: 



oxp I ~^tr0 +tr^V^ 



(16) 



d d 
This is because operator A = V' tttt ^ s more complicated, than the previously considered operator T . 



dtp 1 

Its components do not commute, forming a non-abelian GL(n) algebra 



dtp 1 



A) At - Af A) = 5*Aj - 5\A k 



(17) 



As a result of this, (16) breaks down already in the 2nd order of perturbation theory in t: 



I 



Xi A3 \k 21 , 2i 2k 2j 21 , 2i 2k 21 2j 



f l = 



3(tri 2 ) 2 + (tri) 2 -iVtri 2 



(18) 
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i.e. two additional terms appear at the r.h.s. The leading contribution 3 ^tr A J is what one would expect, 

if identity (16) was true. The other two terms are due to non-abelian nature of operators Aj and can be 
considered as subleading contributions. Following this line of thinking, we state that (16) does not hold 
exactly, but holds approximately and can be considered as a 0-th approximation to a correct identity: 

J d(f> exp ^— -i-tr cj) 2 + corrections + tr 4>A^J = cxp ^ trA 2 ^j (19) 

In the next section we will show, that corrections are suppressed by powers of t. Having this in mind, we 
use (19) to obtain an (approximate) Kontsevich-Hurwitz matrix integral: 



e tW /2 e tr4, = exp 



i t ri 2 -ytri^ e tr ^« J dcj> exp^-ltr^+tr^i)- ytri) e tr ^ 



where the integral is taken over all n x n Hcrmitian matrices with conventional measure. Using 



exp 



(tr Mi) fty) = cxp (^trMV^) /(V) = f(e M ^) VM, /, 



we finally obtain 



> j <lo exp( jj:vo 2 ■ < 2 u(rV) 



(20) 



Emerging matrix model (20) with exponential potential in the presence of background field ip, despite being 
only approximate, is very interesting. Its equation of motion 



_d_ ( 1 



(-ltr0 2 + e-*/ 2 tr (e*V)) = ~4 + ^ = 



or briefly 4>e * = te Nt ^ 2 ip is a transcendental equation on matrix cj>, solved by Lambert function: 
d ty) = y + y 2 + V + % 4 + . . . = £ ^iT, 2/ = te^V 

m—1 

which - naturally - plays a big role in the still underdeveloped and mysterious theory of the Hurwitz- 
Kontsevich model [14]. Thus, despite written in the abelian approximation, matrix integral (20) captures 
correctly a crucially important feature of the Hurwitz-Kontsevich model and deserves further study. 

4.2 Corrections 

Due to non-commutativity of operators A 1 -, there are non- vanishing corrections at the left hand side of 
(19). We introduce them in the form of additional potential 

V(<j>,t) = a(t) + ai(t) tr<p l + a y (i) tvft tr<^ + a ijk (t) trft ti(fP tr<f> k + ... 
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which appears at the left hand side in 



oxp ( -— tT(f) 2 + V(<p,t) +tr^ij =cxp^tri 2 ) (21) 



Potential V(cj), t) should possess a series expansion in non- negative powers of t 

k=0 



because subleading contributions (like those at the right hand side of (18)) are always nested commutators 
of operators A and, therefore, always have lower A-degree than the leading contribution. To cancel these 
terms, potential V(<p,t) must contain higher powers of t. The simplest terms in V can be found by direct 
calculations, similar to (18). A parametrization, which is more convenient for these direct calculations, is 

e v = (3(t) + (3i(t) tr0 4 + /3y(i) tr^ tr<^' + p ijk (t) tr0* tr<^ tr/ + . . . 



where /3i ll ...,i m are another parameters, which can be easily expressed through c*i lt ... > i m and vice versa. 
Computer experiments show, that 0i lt ...,i m = 0, if i\ + . . . + i m is odd. For even i\ + . . . + i mi we have 



Pit) = 1 - YA < n2 !) + TT^ n2 ( n2 - ^ - 82944" 3 ^ 2 ~ ^ + °^ 



M)= ^-^ na -l) + ^nV-l) a +0(*») 



feW ^ ^ 2 - 6) - ^„(„» - l)(5n* - 6) + 0(t ») 
/W*) = ^-^»(»»-l) + 0(,») 
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W) 

042 (*) 
ft*® 

/3 4 n(i) 

/?32l(*) 



181440 
1 



30240 
40 - 7n 2 



483840 
1 



n + 0(t) 
0(t) 

-0(t) 



18144 
1 



+ 0(t) 



69120 
1 

17280 



n + 0(t) 
n + 0{t) 



/?311l(t) 

ftllll(t) 
/?llllll(*) 



5n 3 - 18n 



414720 
1 



+ 0(t) 



17280 

6 - 5n 2 
138240 

1 

27648 
1 



82944 



+ 0(t) 
+ 0(t) 
n + 0(t) 
+ 0(t) 



These terms cancel all subleading contributions up to t 6 , i.e. they make eq. (21) valid up to the order t 6 . 
4.3 Exact matrix integral 

Above results, obtained by direct computer calculations, reveal a nice structure. To see this structure, 
notice that the lowest term (3{t) seems to exponentiate 

Pit) = 1 - y A n{n* 1) + - if g^V - I) 3 + ■ ■ ■ = e~^~^ 



as well as the next two terms: 



A <*> = 24H 1 "^ 2 - 1 ) + TT^ 



n 2 (n 2 -l) 2 + .. 



= .Lnp-nC" 2 -!)*/^ 



24 



-n(n 2 -l)i/24 



24 



This is clearly a hint: coefficients a are simpler, than coefficients /?. Indeed, the simplification happens if 
we take a logarithm of the above series: 

V(4>,t) = log (/?(*)+ ft (t) tr^ + /3 y (t) tr^ tr<^' + ...) = 



-—n(n 2 - 1) + — tr </» 2 tr <b° - — tr e^tr 6 1 - 

24 V ' 24 V V 24 r V 2880 



1 tr 4 tr 6° + — tr </> 3 tr — tr <^> 2 tr </> 2 + 



720 



960 



+ tr^ 6 tr cjP 

181440 



tr^'^tri/) 1 + T^T^-tr 4 tr (j) 2 — — tr </> 3 tr </> 3 + higher order terms 



30240 



12096 



18144 



with an obvious notation tr 0° = n. Somewhat mysteriously, the first item is the central term of the Vi- 
rasoro algebra with c = t/2. As one can see, the potential is simple - much simpler than one could have 
expected. It does not contain terms with more than two traces (at least up to order (j) e ). 
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We conjecture, that higher order terms have the same structure: 

t °° 
V(<f>,t) = --n(n 2 -l) + Yl aytr^tr^ 

(»,j)>(0,0) 



where the sum is taken over all non-negative i,j except 0, 0. We have 

1 1 



a40 ^ 04 = -5760 a31 = ai3 ^14l0 a22 = _ 960 

aeo = a ° 6 = 36^80 a51=ai5 ^"60i80 a42=a24 = 24l92 * 33 = "Islii 



while ajj with odd i + j vanish. Looking at these numbers, it is easy to recognize that 

B l+j . . 
a ij = or , ■ \ —tt ' 1 + 3 = even positive, 

where B 2n are Bernoulli numbers: 

Bo = — , Ba = , Bp. = , B% = , -Bin = — 1 Bin = , . . . , 

6' 30' 42' 8 30' 66' 2730' ' 

generated by 

EB k z k z z ^ B k z k z (z 

— rr~ = 1 + ~ or > — rr~ = x coth - 

kl e z -l 2 ^ k\ 2 V2 

fc=2 fc=2 

Thus our conjecture is that 

jd* exp(-±ta^-±»(» a -l) + g ^yf±f tr 0* tr ^ + tr <f>A 

V i+j>2 



exp(-trA 2 ) (22) 



and exact Hurwitz-Kontsevich matrix integral is 



-I 



d(j) exp 



2t 



-tr< 



h^ 2 - 1) + E 



i+i>2 



(-1) J gi+j 

2(i + j) t!j! 



\ 



tr0 l tr^' + tr (e^- n */ 2 V) 



(23) 



/ 



As one can see, the Gaussian part — tr</> 2 /(2t) is indeed dominating in the small-t limit, because the 
Bernoulli part of potential is of order t° and its constant part is of order t. To emphasize the i°-nature 
of Bernoulli series as a " quasiclassical correction" one can simply move it to the integration measure, by 
summing up the series: 



E 

i,j=0 

i+j>2 



(-1) 3 " B. 
2(* + j) 



tr ft tr (jP = 



= V 



Br, 



' 2m ■ ml 

m=2 



tr (4> <8> I - I <g> (/>) m = -tr log 



sinh 



2 
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and 



exp 



E 

i,3 = 

d+j>2 



B i+j 

2(i + j) 



sinh I 



tr <j) 1 tr (fs? = dct 



1/2 



^ 4>®i-i®4> \^ 



so that 




(24) 



Note that, as usual for GKM [3, 5], this integral does not depend on n if expressed as a function of pk = trip k . 

It may be convenient to shift </> *—> <fr+nt/2. The dcterminantal part of the measure is obviously invariant 
under this shift, so we obtain 



Zhk — 



sinh 



det 



M \ ( 



2 ) 



expl--tr« 



This matrix model is, of course, of eigenvalue type. When expressed through eigenvalues A* of matrix </>, 
the determinant part of the measure takes form 



det 



'sinh( ^V^ ) 



1/2 



n 



(A 4 -A 3 )/2 _ p {\j—\i)/2 



1 — n 



Aj — Xj 



= exp 



E A 0n 

i=i / i<j 



e x% - e x i 
Aj — A, 



and the angular integration can be done with the help of the Itzykson-Zuber formula [4]: 

det b exp (e Xa uu b ) 



[dU] exp (tr (e^t/" 1 )) = =- 

1 Li 



where are eigenvalues of matrix tp. Under the integral sign, it is possible to substitute 



det exp [e Xa uJb) ► exp I ^ e Ai Wi 
\i=i 

In this way we represent Zhk as n-fold integral over eigenvalues: 



z<j J \ i— 1 z— 1 i—1 / 



(25) 



Eq.(24) and its eigenvalue representation (25) is the corrected form of the naive matrix model (20). This 
is an inspiring formula, with many ingredients parallel to intriguing observations about Zhk in [14], but 
further discussion remains beyond the scope of the present paper. Instead we describe in the next subsection 
an alternative representation of Zhk in the form of a discrete matrix model [5]. 
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4.4 Character expansion 



On the space of all polynomials in variables pk, the Hurwitz operator Wq acts in a sophisticated way, 
by mixing different monomials. However, under this sophisticated action certain polynomials - naturally 
called eigenfunctions of Wo _ map into a multiple of themselves. Actually, cigenfunctions \R are characters 
of irreducible representations of GL(n), labeled by their signatures 

R= (Ai,A 2 ,...A m ), with Ai > A 2 > ... > A m 

also known as partitions or Young diagrams. In this section, we will use two classical formulas from the 
theory of characters. 

The first formula defines the eigenvalues of the operator Wo: 

m 

W oX r = C rX r, C R = J2 x i ( x i - 2i + 1) (26) 
We give a proof of this eigenvalue formula in the Appendix in s.4.5. The first several eigenfunctions are 



Y(71 = 1 




Xi =Pi, 


W oX i = 


1 1 2 

X2 = 2^2 + 2P1, 


W0X2 = +2X2 


1 1 2 

Xi,i = -^P2 + -jVi, 


WoXi,i = -2xi,i 


11 1 , 

X3 = gP3 + -jP2Pl + gP?, 


WoX3 = 6x3 


1 1 , 

X2,l = -gP3 + gPl, 


WoX2,i = 


11 1 , 

Xl,l,l = gP3 - 2"P2Pl + gPl, 


WoXi.1,1 = -6X1,1,1 


1 1 2 1 1 2 1 4 

X4 = ^4 + gP 2 + gPlP3 + ^P2Pt + ^Pl, 


W oX 4 = 12X4 


1 1 , 1 2 1 . 

X3,l = -^P4 - gP 2 + ^P2Pi + gPl, 


^0X3,1 = 4x3,1 


1 2 1 1 I 
X2,2 = ^2 - gPlPS + jgPT, 


WoX2,2 = 


1 1 2 1 2 1 4 

X2,l,l - ^P4 - gP2 - -^P2P{ + gPT, 


W X2,1,1 = -4^2,1,1 


1 1 1 1 2 1 

Xi,i,i,i = -4P4 + gP 2 + 3P1P3 - -^>2P{ + ^Pi, 


^0X1,1,1,1 = -12xi,i,i,i 



The second formula states, that e Pl is decomposed into eigenfunctions with the following coefficients: 
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II (Xi-Xj-i + j) 

e pi = E d ^R, d R = (27) 

R Y[{\ t + m-i)\ 

1=1 

It is well known as the hook formula [22] . As we show in s.4.5, it is a particular case of more general Cauchy 
identity (30). Up to the simple n-dependent factors d R is the dimension of representation R of GL(n): 

1=1 v ' 

The first several coefficients dn and corresponding dimensions dim# are 



d~ — 1 


r\ i m ™ — ' 

\11L1L0 — J_ 


di = 1, 


dimi = n 


d 2 = \, 


n(n + l) 
dim 2 = 


di,i = \, 


rc(n-l) 
dimi,! = 


d3 = l' 


n(n + l)(n + 2) 

dim 3 = 

6 


<h,i = \, 


n{n 2 - 1) 
dim 2i i = 




dim,,,,^" 1 "- 1 ''"- 21 


rf4= 2l' 


n(n+ l)(n + 2)(n + 3) 
4 = 24 




dim 31 = n(n + 2)(n2 - 1) 


d2,2 = ^ 


n 2 (n 2 - 1) 
dim 2:2 = ^ 


d2,i,i = \, 




di,M,i = ^. 


dim 11 , 1 -" ( "- 1)( "- 2)( "- 3) 



Eqs. (26) and (27) imply the following character expansion for Zhk ■ 
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Zhk = eWR ' 2 dR * R 

R 



This sum over partitions, i.e. a sum over representations of GL(n), can be regarded as a discrete matrix 
model [5] and can be handled by the methods of [23] . Its equivalence to the matrix model (24) is somewhat 
obscure and will be discussed elsewhere. A few first terms of the character expansion (28) are 



(28) 



Z H K(p\t) 



l + Xi + ^e*X2 + ^e *xn + ie 3t X3 + ^X2i + 3 *xm + 
^e 6t X 4 + ^e 2t X 3i + ^X22 + ^e- 2t X2ii + ^e- 6t xmi + . . . 



or directly through the time-variables 

ZHI<(p\t) = 



1 _ L, i \ ,,2 / J_, .",/ _ _L, 1 „,„, j. [ _L„6t _ J__ 2t , J L p -2i i J_ p -6t 1 ,,2 



1+1 +\ n ( 1 q+ 1 



4 4 / ^ V 12 12 I y y V 192 64 48 64 192 ' F2 



—e 6t - — + —e- 6t ) p 3Pl + ( — e 3t + - + — e" 3 * ) p 3 + ( — e 6 * + — e 2t - — e" 2 * - — 
72 36 72 / ^ V 36 9 36 / Pl \ 96 32 32 96 



f ;[ 6t _i_ JL„2t j 1 [_ JL c -2t j (>l } , 



-e m + —e zz + + — e- zl + e~ m + . . 

V576 64 144 64 576 ' 

Note, that in this way we obtain a series in p-variables, but summed in all orders of the i-variablc. 

4.5 Appendix: GL(n) characters as eigenfunctions of W 

In this review section our starting point is the free-fermion Wick theorem [5, 24], 



m m 



( n n ^(w) ) = det < faMvi) ) 



i=l i=l 



- a generalized form of the Fay identity, KP-equations and Shottky relations,- which in the flat coordinates 
on the Riemann sphere reduces to the elementary Cauchy identity: 

A(x)A(y)l[^— = det— ^— , 
where A{x) = Ili<j( x i ~ x j)- Equivalently, it can be written as 

TT = a/ hn ^ det — ( 29 ) 

The series expansion of the right hand side has the form 

1,1 _ 
4f i =1^XRXR 



A(x)A(y) i,j 1 - XiVj 



R 
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where sum goes over irreducible representations R = (Ai > A2 > . . . A TO ) of GL(n), and 

Het ( r ^i+m-i\ Apf l n A 4 +m-»\ 

X-R = XAi,...,A m = ' Xfl = XAi,...,A m = / m -i\ 

are the characters. Since characters are symmetric functions in Xi and yj, by inverse Miwa transform they 
can be expressed through the time variables ak — x\ + . . . + xj^ and bk = y k + . . . + (as in [14] we use 
the time-variables which differ from conventional ones, tk, by a fc-factor: pk = ktk). Expressed through a& 
and bk, (29) takes the form 

exp ( k akbk I = ^ XA 1 ,...,A m (a)xA 1 ,...,A m (&) = ^2xR(a)XR(b) (30) 

\fe=l / Ai>A 2 >...>A m R 

In this form, it is also known as Cauchy identity. In particular case of ak = Sk,i, we recover (27): 

exp(6i) = ^Xfl(Mxfl(6) (31) 

R 

where coefficients du are given by characters 

m 

n (Ai-Aj-i+j) 

eta = Xfl(^fe.i) = 

n (Ai+m-i)! 
»=i 

The last equality requires a straightforward algebraic verification, which we do not include in this paper. 

Apart from the two forms of Cauchy identity, a lot of equally explicit formulas are known for characters, 
see, for example, [3, 4, 22]. However, Cauchy identity is quite enough for our purpose. If we put a k = 
x\ + . . ■ + x k n and bk = Pk, then (30) turns into a "generating function" for m-index characters w.r.t to the 
variables x\, . . . , x m : 

(™x k + ...+x k m \ x . det mxm {xf +m ~ i ) 

\fe=l K / A 1 >A 2 >...>A m clet mxm ^- j 

The simplest generating function corresponds to m = 1 

X{x) = exp I E — p k J = 2J ^ fe Xfc 

\fe=l / k 

and the next-to-simplest corresponds to m = 2: 

X (x, y) = exp I £ ^— Pk = ^ — j Xk,i 

\k=l I k>l v y 7 
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Using these generating functions as a definition of polynomials XAi,...,A m5 we will now prove, that they are 
indeed cigcnfunctions of the Hurwitz operator. Notice, that 



d 



o X{ x l> • • • ! x m) — ^ ] XO^li • • • i 



i=l 



and 



— x(xi, ■ ■ -,x m ) = ^Paxl 1 x(xi, ■■ -,x m ) 



0=1 



Therefore 



00 m 



W X(X1, ■ ■ ■ ,X m ) = £ Y,PaPbX a i +h + £ E ft+i) 1 " 1 ' X(^l,---^m) 
\a,b=li=l a,6=li,j = l / 



00 m 00 m 00 x — X' 1 ■ 

£ £ PaP6< +& +££(*" l) Ps < + E E Pa ' ' ' ^ 

a,fc=li=l s=2i=l s=2i=£j x i x j 



X(.Xl , ■ ■ • , &ra) 



\ ^^"i ^"^3 



d d 



X{xi 1 ■ ■ ■ iXn 



We have just proved, that generating functions satisfy 



W x(xi, ■ ■ -,x m ) 



dx 2 



d 

dxi 



x{xi , ■ • ■ , x m ) 



(33) 



The operator on the right hand side is a Hamiltonian of Calogcro-type dynamical system. The first term in 
square brackets is the kinetic energy, while the second term represents interaction. Instead of diagonalizing 
Wo, we can do it for this Hamiltonian, which is much simpler. Let us show, that 



x 2 d 2 \ x;.v j 



i=i 



d_ _d_ 

1 dx 2 ' ^ Xi — Xj V dxi dx 

in 



2^ X i f)^2 



Het („Ai+m-i\ j + , 

^ ct mxm y-^j J ^ ^ clJ mxm y^j 

ueu mX m y^j ) 



\i-\-m—i\ 



det mX m (Xj ) 



where C\ 1> ... t \ m is a number (does not depend on x). This is obvious for m = 1 

x 2 ^ = (fc 2 -^ 



and almost obvious for to = 2: 



2 9 2 , 2 a 2 



2xy ( d d 



dx 2 dy 2 x — y \ dx dy 



x k+1 y l - x l y k+1 \ = ( k 2 + l 2 _ k _ 3l jf xk+1 v l - xl v k+ ' 



x-y 



x-y 



For higher to, the proof is a straightforward algebraic exercise, which we present here in full detail. To 
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begin with, wc show that both A = det (x™ l ) and A = det [x- i+m ') are eigenfunctions of kinetic energy. 
Indeed, making use of explicit formulas for both determinants 



A = det (x 



X{ -\-m—i 
3 



)= E(-!) H - 



\a\ ^"1 +m — <ri x ^<r m +m—cr m 



and 



A= det (xf- 1 ) 



E (-i) kl < 



we have 



E = E w E + - - -0 + m - ^ - 1) 4 

i=l 1 cr£S m i=l 



• • • 



The sum over i, which appears in the right hand side, is invariant under permutations and factors out: 

,9 2 A 



E^T^T = + m - + m - i - l) A 



Oxf 



Similarly 



E 



,d 2 A 



1 dx 2 
l a * aes. 



E (~ 1 ) H E ( m ~ ( m ~ Ui - l ) x T^ ai ■ ■ ■ x 2~° m = E ( m ~ ( m ~ 1 ~ l ) A 

i=l 



i=l 



So A and A are eigenfunctions of kinetic energy. By straightforward differentiation, we obtain 



E 



d 2 f A\ 1 A 9 d 2 A A^ 9 <9 2 A 2 ^ 9 <9A5A 2 ^ 9 ~ <9A 5 A 
" x 



^ 1 dx 2 \ A ) 1 dx 2 A 2 f-f 4 9a;, 2 A 2 f-' 1 ftr, A 3 f-' 1 5x l 5a; 

2=1 1 \ / 2 — 1 1 2—1 * 2 — 1 2—1 



-+A3E^ 



= E(A i + 2m-2*-l) ---J^J 

i—\ i—i i—i 

Including interactions and considering the full Hamiltonian, we obtain 



d 2 



j2 x 2 — h XiXj [ — — 

2=1 ^7^3 ~~ "^3 ^^3 



d d 



Am A 

-r-EA 4 (A l + 2m-2z-l) - 
^ i=l ^ 



2xjXj / j_ 9A _ _A_ 9A \ _ _2_ ^ ^ <9A <9A + J_y x 2 7 L dA dA 

ijtj Xi -Xj I A dxi A 2 dxi J A 2 j=1 1 fej da^ A 3 i=1 1 fei ftcj 



The sum over j is easily evaluated 
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E 



LJ- = Xi -2m + 2 + — 



and, after cancelation of terms, we obtain 



m 2 <9 2 

S X i "^~2 + E 



-( d 




■ \dxi 


dxjj 



A m A 
-r-EA i (A J + 2m-2«-l) - 
^ i=l A 



-X,*, ^-2m + 2 + — — j _ + _ X^A— _ 

m / 1 9A A f ~ \ A 

= g a, (-2m + 2) ^_____j = (-2m + 2) (deg A - deg a) - 

_ m m 

By definition, deg A = E (^» + * — 1) an d deg A = E ~ !)■ So, we have proved the identity 



d 2 ^ x i x j f ® ® 
doC' Xi x i V dxi dx^ 



E 



A , .A 

i=i 



(34) 



thus diagonalizing this Hamiltonian. At the same time, via Cauchy identity, we have proved the dual result 



^ ^( a + & )P«^^^+a^a+6^^JxA 1 ,...,A m (p) =5^Ai(Ai-2i+l) XA 1 ,...,A m (p) (■•-!•".) 



(3) 

thus diagonalizing the Wq operator with quadratic eigenvalues 



C Xl ,...,x m =Y, x i( x i- 2i + 1 ) ( 36 ) 

i=l 

Above calculation is somewhat tedious, but important: it shows, that the Wq "cut-and-join" operator, 
defining the generating function of Hurwitz numbers, can be rewritten as a Calogero-type Hamiltonian of 
multi-particle system in Id with coordinates Xi and pairwise interactions between particles. 



5 Conclusion 

In this paper we considered three formulas, relating the two a priori unrelated kinds of objects: partition 
functions of the matrix models and generators of the W^-algebra. 

Eq.(4) provides a one more representation for Hcrmitian matrix model - the main personage of all 
matrix-model studies. 

Eq.(3) describes a mush less investigated version of the same model, with extra background field tp. 
The two VF-operators in the two formulas (3) and (4) act on the two conjugate sets of time-variables in 
Hcrmitian model. It would be interesting to extend these expressions to non-Gaussian, say, DV phases of 
the theory. 

Finally, eqs.(13) and (24) provide an inspiring matrix-model representation for the Hurwitz-Kontsevich 
tau-function - a recently discovered new link between combinatorics and integrability theory. It would be 
very interesting to explicitly describe the intriguingly sophisticated Virasoro constraints for this partition 
function as the Ward identities for this matrix model. 



22 



Acknowledgements 



We are indebted to A.Mironov for important comments. Our work is partly supported by Russian Federal 
Nuclear Energy Agency and the Russian President's Grant of Support for the Scientific Schools NSh- 
3035.2008.2, by RFBR grant 07-02-00645, by the joint grants 09-01-92440-CE, 09-02-91005- ANF, 09-02- 
93105-CNRS and by the NWO project 047.011.2004.026. The work of Sh.Shakirov is also supported in 
part by the Moebius Contest Foundation for Young Scientists and by the Dynasty Foundation. 

References 

[1] E.Brezin, C.Itzykson, G.Parisi and J.-B.Zuber, Planar diagrams, Comm. Math. Phys. 59 (1978) 35; 
D.Bcssis, A new method in the combinatorics of the topological expansion, Comm. Math. Phys. 69 (1979) 
147; 

D.Bcssis, C.Itzykson and J.-B.Zuber, Adv. Appl. Math. 1 (1980) 109; 

M.-L. Mehta, A method of integration over matrix variables, Comm. Math. Phys. 79 (1981) 327; Ran- 
dom Matrices, 2nd edition, Acad. Press., N.Y., 1991; 

D.Bcssis, C.Itzykson and J.-B.Zuber, Quantum field theory techniques in graphical enumeration, 
Adv.Appl.Math. 1 (1980) 109; 

A.Migdal, Loop equations and 1/N expansion, Phys. Rep. 102 (1983) 199; 

F.David, A Model of Random Surfaces with Nontrivial Critical Behavior, Nucl. Phys. B257 [FS14] 

(1985) 45, 543; 

J. Ambjorn, B. Durhuus and J. Frohlich, Diseases of Triangulated Random Surface Models, and Pos- 
sible Cures, Nucl. Phys. B257 [FS14] (1985) 433; 

V. A. Kazakov, I. K. Rostov and A. A. Migdal, Critical Properties of Randomly Triangulated Planar 
Random Surfaces, Phys. Lett. 157B (1985) 295; 

D. Boulatov, V. A. Razakov, I. R. Rostov and A. A. Migdal, Possible Types Of Critical Behavior And 
The Mean Size Of Dynamically Triangulated Random Surfaces, Phys. Lett. B174 (1986) 87; Analytical 
and Numerical Study of the Model of Dynamically Triangulated Random Surfaces, Nucl. Phys. B275 
[FS17] (1986) 641; 

V. Razakov, Ising model on dynamical planar random lattice : exact solution, Phys. Lett. A 119 

(1986) 140, The appearance of matter fields from quantum fluctuations of 2D -gravity, Mod. Phys. Lett. 
A4 (1989) 2125; 

E. Brezin and V. Razakov, Exactly Solvable Field Theories Of Closed Strings, Phys. Lett. B236 (1990) 
144; 

D. Gross and A. Migdal, Nonperturbative Solution of the Ising Model on a Random Surface, Phys. Rev. 
Lett. 64 (1990) 127; A Nonperturbative Treatment Of Two- Dimensional Quantum Gravity, Nucl. Phys. 
B340 (1990) 333; 

A.Mironov and A.Morozov, Phys.Lett. B252(1990) 47-52; 

F. David, Loop Equations and Nonperturbative Effects In Two- Dimensional Quantum Gravity, 
Mod.Phys.Lett. A5 (1990) 1019; 

J.Ambjorn and Yu.Makeenko, Mod.Phys.Lett. A5 (1990) 1753; 

L. Alvarez-Gaume, Random surfaces, statistical mechanics, and string theory, Lausanne lectures, 1990; 
A. Levin and A.Morozov, On the Foundations of the Random Lattice Approach to Quantum Gravity, 
Phys.Lett. 243B (1990) 207-214; 

H.Itoyama and Y.Matsuo, Phys.Lett. B255 (1991) 202; 
P.Ginsparg, Matrix Models of 2d Gravity, hcp-th/9112013; 

A.Marshakov, A.Mironov and A.Morozov, On the Equivalence of Topological and Quantum 2d Gravity, 
Phys.Lett. 274B (1992) 280-288, hcp-th/9201011; 



23 



V.Kazakov and A.Migdal, Induced QCD at large N, Nucl.Phys. B397 (1993) 214-238,1993, hep- 
th/9206015; 

P. Di Francesco and C. Itzykson, A Generating Function for Fatgraphs, Annales Poincare Phys.Theor. 
59 (1993) 117-140, hcp-th/9212108; 

J.-M. Daul, V.A. Kazakovand I.K. Kostov, Rational theories of 2D gravity from the two-matrix model, 
Nucl. Phys. B409 (1993) 311; 

M. Staudacher, Combinatorial solution of the two matrix model, Phys. Lett. B305 (1993) 332, hep- 
th/9301038; 

J.Ambjorn, L. Chekhov, C.F.Kristjansen and Yu.Makeenko, Matrix Model Calculations beyond the 
Spherical Limit, Nucl.Phys. B404(1993) 127-172, Erratum B449 (1995) 681, hep-th/9302014; 
P.Di Francesco, P. Ginsparg and J. Zinn- Justin, 2D Gravity and Random Matrices, Phys. Rep. 254 
(1995) 1-133, hcp-th/9306153; 

S.Kharchev, A.Marshakov, A.Mironov and A.Morozov, Generalized Kazakov-Migdal-Kontsevich 
Model: group theory aspects, Int. J. Mod. Phys. A10 (1995) 2015, hep-th/9312210; 
B.Eynard, Large Random Matrices: Eigenvalue Distribution, hcp-th/9401165; 

V.Kazakov, M. Staudacher and Th.Wynter, Character expansion methods for matrix models of dually 
weighted graphs, Commun.Math.Phys. 177 (1996) 451-468, hcp-th/9502132; 

M.Adler, A.Morozov, T.Shiota and P. van Moerbeke, New matrix model solutions to the Kac-Schwarz 
problem, Nucl.Phys.Proc.Suppl. 49 (1996) 201-212, hcp-th/9603066; 

G. Akemann, Higher genus correlators for the Hermitian matrix model with multiple cuts, Nucl.Phys. 
B482 (1996) 403-430, hcp-th/9606004; 

G. Akemann, P.H.Damgaard, U.Magnca and S.Nishigaki, Universality of random matrices in the mi- 
croscopic limit and the Dirac operator spectrum, Nucl.Phys. B487 (1997) 721-738, hep-th/9609174; 
Multicritical microscopic spectral correlators of Hermitian and complex matrices, Nucl.Phys. B519 
(1998) 682-714, hcp-th/9712006; 

A. Marshakov, M.Martellini and A.Morozov, Insights and puzzles from branes: J^-D SUSY Yang-Mills 
from 6-D models, Phys.Lctt. B418 (1998) 294-302, hcp-th/9706050; 

T.Guhr, A.Mueller-Groeling and H.A.Weidenmueller, Random Matrix Theories in Quantum Physics: 
Common Concepts, Phys. Rep. 299 (1998) 189-425, cond-mat/9707301; 

V. A. Kazakov, I. K. Kostov and N. A. Nekrasov, D-particles, matrix integrals and KP hierarchy, Nucl. 
Phys. B 557 (1999) 413, hcp-th/9810035; 

B. Eynard, Random Matrices (2000), http://www-spht.cea.fr/articles_k2/t01/014/publi.pdf; 

P. Forrester, N. Snaith and J. Verbaarschot, Developments in Random Matrix Theory, J. Phys. A36 
2859-3645, cond-mat/0303207; 

S.Alexandrov, V.Kazakov and D.Kutasov, Nonperturbative effects in matrix models and D-branes, 
JHEP 0309 (2003) 057, hcp-th/0306177; 

V.Kazakov and I. Kostov, Instantons in noncritical strings from the two matrix model, hep-th/0403152; 
G. Akemann, Y.V.Fyodorov and G. Vernizzi, On matrix model partition functions for QCD with chem- 
ical potential, Nucl.Phys. B694 (2004) 59-98, hcp-th/0404063; 

P.Di Francesco, 2D Quantum Gravity, Matrix Models and Graph Combinatorics, math-ph/0406013; 
A.Morozov, Challenges of matrix models, hep-th/0502010; 

J.Harnad, A.Orlov, Fermionic construction of tau functions and random processes, Physica D235 
(2007) 168-206, arXiv:0704.1157 

[2] A. Gerasimov, A. Marshakov, A. Mironov, A. Morozov, and A. Orlov, Matrix Models of Two- 
Dimensional Gravity and Toda Theory, Nucl. Phys. B357 (1991) 565618; 

Yu. Makeenko, A. Marshakov, A. Mironov and A. Morozov, Continuum versus discrete Virasoro in 
one-matrix models, Nucl.Phys. B356 (1991) 574; 



24 



S.Kharchev, A.Marshakov, A.Mironov and A.Morozov, Generalized Kontsevich Model Versus Toda 
Hierarchy and Discrete Matrix Models, Nucl.Phys. B397 (1993) 339-378, hep-th/9203043; 
S.Kharchev, A.Marshakov, A.Mironov, A.Morozov and S.Pakuliak, Conformal Matrix Models as an 
Alternative to Conventional Multi-Matrix Models, Nucl.Phys. B404 (1993) 717-750, hcp-th/9208044 

[3] M. Kontsevich, Intersection theory on the moduli space of curves, Funkts. Anal. Prilozh., 25:2 (1991) 
5057; 

M. Kontsevich, Intersection theory on the moduli space of curves and the Airy function, 
Comm.Math.Phys. 147 (1992) 1-23; 

S.Kharchev, A.Marshakov, A.Mironov, A.Morozov and A.Zabrodin, Unification of All String Models 
with c < 1, Phys. Lett. B275 (1992) 311-314, hcp-th/9111037; Towards unified theory of 2d gravity, 
Nucl.Phys. B380 (1992) 181-240, hcp-th/9201013; 
M.Adlcr and P. van Moerbeke, Comm.Math.Phys. 147 (1992) 25; 

P.Di Francesco, C.Itzykson and J.-B.Zuber, Polynomial averages in the Kontsevich model, 
Comm.Math.Phys. 151 (1993) 193-219, hep-th/9206090 

[4] Harish-Chandra, Am. J. Math. 79 (1957) 87; 

C.Itzykson and J.-B.Zuber, J.Math.Phys. 21 (1980) 411; 
J.Duistermaat and G.Heckman, Invent. Math. 69 (1982) 259; 

A. Hietamaki, A.Morozov, A.Niemi and K.Palo, Geometry of N= 1/2 supersymmetry and the Atiyah- 
Singer index theorem, Phys.Lctt. B263 (1991) 417-424; 

Super symplectic geometry of super symmetric quantum field theories, Nucl.Phys. B377 (1992) 295-338; 
M.Bowick, A.Morozov and D.Shevitz, Reduced unitary matrix models and the hierarchy of tau func- 
tions, Nucl.Phys. B354 (1991) 496-530; 

I.Kogan, A.Morozov, G.Semenoff and H.Weiss, Area Law and Continuum Limit in "Induced QCD", 
Nucl.Phys. B395 (1993) 547-580, hep-th/9208012; Continuum Limits of "Induced QCD": Lessons of 
the Gaussian Model at d=l and Beyond, Int.J.Mod.Phys. A8 (1993) 1411-1436, hep-th/9208054; 
A.Morozov, Pair Correlator in the Itzykson-Zuber Integral, Mod. Phys. Lett. A7 (1992) 3503-3508, 
hep-th/9209074; 

S.Shatashvili, Correlation Functions in The Itzykson-Zuber Model, Comm.Math.Phys. 154 (1993) 421- 
432, hcp-th/9209083; 

A. Mironov, A.Morozov and G.Semenoff, Unitary matrix integrals in the framework of General- 
ized Kontsevich Model. I. Brezin- Gross- Witten Model, Int.J.Mod.Phys. All (1996) 5031-5080, hep- 
th/9404005; 

B. Eynard, A short note about Morozov's formula, math-ph/0406063; 

A. Ferrer, B. Eynard, P. Di Francesco and J.-B. Zuber, Correlation Functions of Harish-Chandra 
Integrals over the Orthogonal and the Symplectic Groups, J. Stat. Phys. 129 (2009) 885-935, math- 
ph/0610049; 

M.Bergere and B. Eynard, Some properties of angular integrals, arXiv:0805.4482 

[5] A.Morozov, Integrability and Matrix Models, Phys.Usp. 37 (1994) 1-55, hep-th/9303139; Matrix Models 
as Integrable Systems, hep-th/9502091; 

A.Mironov, Matrix Models vs. Matrix Integrals, Theor.Math.Phys. 146 (2006) 63-72, hep-th/0506158 

[6] R.Dijkgraaf and C.Vafa, Matrix Models, Topological Strings, and Supersymmetric Gauge Theories, 
Nucl.Phys. B644 (2002) 3-20, hep-th/0206255; On Geometry and Matrix Models, Nucl.Phys. B644 

(2002) 21-39, hcp-th/0207106; A Perturbative Window into Non-Perturbative Physics, hep-th/0208048; 
L.Chekhov and A.Mironov, Matrix models vs. Seiberg- Witten/ Whitham theories, Phys. Lett. B552 

(2003) 293-302, hep-th/0209085; 

R.Dijkgraaf, S.Gukov, V.Kazakov and C.Vafa, Analysis of Gauged Matrix Models, Phys. Rev. D68 



25 



(2003) 045007, hcp-th/0210238; 

V.Kazakov and A.Marshakov, Complex Curve of the Two Matrix Model and its Tau-function, J.Phys. 
A36 (2003) 3107-3136, hcp-th/0211236; 

H. Itoyama and A.Morozov, The Dijkgraaf-Vafa prepotential in the context of general Seiberg-Witten 
theory, Nucl.Phys.B657 (2003) 53-78, hcp-th/0211245; Experiments with the WDVV equations for 
the gluino- condensate prepotential: the cubic (two-cut) case, Phys.Lett. B555 (2003) 287-295, hep- 
th/0211259; Calculating Gluino- Condensate Prepotential, Prog.Theor.Phys. 109 (2003) 433-463, 
hep-th/0212032; Gluino- Condensate (CIV-DV) Prepotential from its Whitham-Time Derivatives, 
Int.J.Mod.Phys. A18 (2003) 5889-5906, hcp-th/0301136; 

S.Naculich, H.Schnitzer and N. Wyllard, Matrix model approach to the N~2 U(N) gauge theory with 
matter in the fundamental representation, JHEP 0301 (2003) 015, hep-th/0211254; 
B.Feng, Geometric Dual and Matrix Theory for SO/Sp Gauge Theories, Nucl.Phys. B661 (2003) 
113-138, hep-th/0212010; 

I. Bena, S.de Haro and R.Roiban, Generalized Yukawa couplings and Matrix Models, Nucl.Phys. B664 
(2003) 45-58, hcp-th/0212083; 

Ch.Ann, Super symmetric SO(N)/Sp(N) Gauge Theory from Matrix ModeV.Exact Mesonic Vacua, 
Phys.Lett. B560 (2003) 116-127, hcp-th/0301011; 

L.Chekhov, A.Marshakov, A.Mironov and D.Vasiliev, DV and WDVV, hcp-th/0301071; Complex Ge- 
ometry of Matrix Models, Proc. Stcklov Inst. Math. 251 (2005) 254, hcp-th/0506075; 
A. Dymarsky and V. Pcstun, On the property of Cachazo-Intriligator-Vafa prepotential at the ex- 
ternum of the superpotential, Phys.Rcv. D67 (2003) 125001, hcp-th/0301135; 

Yu.Ookouchi and Yo.Watabiki, Effective Superpotentials for SO/Sp with Flavor from Matrix Models, 
Mod.Phys.Lctt. A18 (2003) 1113-1126, hcp-th/0301226; 

H.Itoyama and H.Kanno, Super eigenvalue Model and Dijkgraaf-Vafa Proposal, Phys.Lett. B573 (2003) 
227-234, hep-th/0304184; Whitham Prepotential and Superpotential, Nucl.Phys. B686 (2004) 155-164, 
hcp-th/0312306; 

M.Matone and L.Mazzucato, Branched Matrix Models and the Scales of Super symmetric Gauge The- 
ories, JHEP 0307 (2003) 015, hep-th/0305225; 

R.Argurio, G.Ferretti and R.Heise, An Introduction to Supersymmetric Gauge Theories and Matrix 
Models, Int.J.Mod.Phys. A19 (2004) 2015-2078, hep-th/0311066; 

M.Gomez-Reino, Exact Superpotentials, Theories with Flavor and Confining Vacua, JHEP 0406 (2004) 
051, hep-th/0405242; 

K.Fujiwara, H.Itoyama and M.Sakaguchi, Supersymmetric U(N) Gauge Model and Partial Breaking of 
N=2 Supersymmetry, Prog.Theor.Phys. 113 (2005) 429-455, hep-th/0409060; Partial Breaking of N=2 
Supersymmetry and of Gauge Symmetry in the U(N) Gauge Model, Nucl.Phys. B723 (2005) 33-52, 
hep-th/0503113; Supersymmetric U(N) Gauge Model and Partial Breaking of N=2 Supersymmetry, 
Prog.Theor.Phys. Suppl. 164 (2007) 125-137, hcp-th/0602267; 

Sh.Aoyama, The Disc Amplitude of the Dijkgraaf-Vafa Theory: 1/N Expansion vs Complex Curve 
Analysis, JHEP 0510 (2005) 032, hcp-th/0504162; 

D.Bcrenstcin and S.Pinansky, Counting conifolds and Dijkgraaf-Vafa matrix models for three matrices, 
hep-th/0602294 

[7] T.J.Hollowood, A.Iqbal and C.Vafa, Matrix Models, Geometric Engineering and Elliptic Genera, JHEP 
0803 (2008) 069, hcp-th/0310272 

[8] A.Zabrodin, New applications of non-hermitian random matrices, cond-mat/02 10331; 

P. Wiegmann and A. Zabrodin, Large N expansion for normal and complex matrix ensembles, hep- 
th/0309253 ; 



26 



R. Teodorescu, E. Bettelhcim, O. Agam, A. Zabrodin and P. Wiegmann, Normal random matrix 
ensemble as a growth problem, Nucl.Phys. B704 (2005) 407-444, hcp-th/0401165 

[9] A.Marshakov and N.Nekrasov, Extended Seiberg-Witten Theory and Integrable Hierarchy, JHEP 0701 
(2007) 104, hep-th/0612019 

[10] A.Alexandrov, A.Mironov and A.Morozov, Partition functions of matrix models as the first special 
functions of string theory. I: Finite size Hermitean 1-matrix model, Int.J.Mod.Phys. A19 (2004) 4127, 
hep-th/0310113; Unified description of correlators in non- Gaussian phases of Hermitean matrix model, 
Int.J.Mod.Phys. A21 (2006) 2481-2518, hep-th/0412099; Solving Virasoro Constraints in Matrix Mod- 
els, Fortsch.Phys. 53 (2005) 512-521, hcp-th/0412205; 

A. Alexandrov, A.Mironov, A.Morozov and P.Putrov, Partition Functions of Matrix Models as the 
First Special Functions of String Theory. II. Kontsevich Model, arXiv:0811.2825 

[11] B.Eynard, All genus correlation functions for the hermitian 1-matrix model, JHEP 0411 (2004) 031, 
hep-th/0407261; Large N expansion of the 2-matrix model, JHEP 0301 (2003) 051, hcp-th/0210047; 
Master loop equations, free energy and correlations for the chain of matrices, JHEP 0311 (2003) 018, 
hep-th/0309036; 

B. Eynard and N.Orantin, Topological expansion of the 2-matrix model correlation functions: diagram- 
matic rules for a residue formula, JHEP 0612 (2006) 026, math-ph/0504058; 

L.Chekhov and B.Eynard, Hermitean matrix model free energy: Feynman graph technique for all gen- 
era, JHEP 0603 (2006) 014, hep-th/0504116; Matrix eigenvalue model: Feynman graph technique for 
all genera, JHEP 0612 (2006) 026, math-ph/0604014; 

B.Eynard, M.Marino and N.Orantin, Holomorphic anomaly and matrix models, JHEP 0706 (2007) 
058, hep-th/0702110; 

N.Orantin, Gaussian matrix model in an external field and non-intersecting Brownian motions, 
arXiv:0803.0705 

[12] J.Schwarz, Superconformal Chern-Simons Theories, JHEP 0411 (2004) 078, hcp-th/0411077; 

A. Basu and J.A.Harvey, The M2-M5 Brane System and a Generalized Nahm's Equation, Nucl.Phys. 
B713 (2005) 136-150, hep-th/0412310; 

J.Bagger and N.Lambert, Modelling Multiple M2's, Phys.Rev. D75 (2007) 045020, hcp-th/0611108; 
Gauge Symmetry and Supersymmetry of Multiple M2-Branes, Phys.Rev. D77 (2008) 065008, arXiv: 
0711.0955; Comments on Multiple M2-Branes, JHEP 0802 (2008) 105, arXiv: 0712.3738; Three- 
Algebras and N=6 Chern-Simons Gauge Theories, Phys.Rev. D79 (2009) 025002, arXiv:0807.0163; 
A.Gustavsson, Algebraic Structures on Parallel M2-branes, arXiv: 0709.1260; Self dual Strings and 
Loop Space Nahm Equations, arXiv: 0802.3456; 
S.Mukhi and C.Papageorgakis, M2 to D2, arXiv: 0803.3218; 

M.Bandres, A.Lipstcin and J.Schwarz, N=8 Superconformal Chern-Simons Theories, arXiv: 
0803.3242; 

A.Morozov, On the Problem of Multiple M2 Branes, JHEP 0805 (2008) 076, arXiv:0804.0913; 
J.Gomis, G.Milancsi, and J.G.Russo, Bagger-Lambert Theory for General Lie Algebras, 
arXiv:0805.1012 v2; 

S.Bcnvenuti, D.Rodriguez-Gomez, E.Tonni and H.Verlinde, N=8 superconformal gauge theories and 
M2 branes, arXiv:0805.1087; 

P.-M.Ho, Y.Imamura and Y.Matsuo, M2 to D2 revisited, arXiv:0805.1202; 

A.Morozov, From Simplified BLG Action to the First- Quantized M-Theory, JETP Lett. 87 (2008) 
659-662, arXiv: 0805. 1703; 

O.Aharony, O.Bergman, D.L.Jafferis and J.Maldacena, N—6 superconformal Chern- Simons-matter 
theories, M2-branes and their gravity duals, JHEP 0810 (2008) 091, arXiv:0806.1218; 



27 



I.A.Bandos and P.K.Townsend, Light-cone M5 and multiple M2-branes, Class. Quant. Grav. 25 (2008) 
245003, arXiv:0806.4777; SDiff Gauge Theory and the M2 Condensate, JHEP 0902 (2009) 013, 
arXiv:0808.1583; 

J.Gomis, D.Rodriguez-Gomez, M. Van Raamsdonk and H.Verlindc, A Massive Study of M2-brane 
Proposals, JHEP 0809 (2008) 113, arXiv:0807.1074; 

J.A.Minahan, W.Schulgin and K.Zarembo, Two loop integr ability for Chern-Simons theories with N=6 
super symmetry, arXiv:0901.1142 

[13] A.Givental, Semisimple Frobenius structures at higher genus, math. AG/0008067; 

A.Alexandrov, A.Mironov and A.Morozov, Instantons and Merons in Matrix Models, Physica D235 

(2007) 126-167, hcp-th/0608228; M-Theory of Matrix Models, hep-th/0605171; 

N.Orantin, Symplectic invariants, Virasoro constraints and Givental decomposition, arXiv:0808.0635; 

[14] A.Hurwitz, Uber Riemann'sche Flachen mit gegebenen Verzweigungpunkten, Math. Ann. 39 (1891) 1- 
61; Uber die Anzal der Riemann'sche Flachen mit gegebenen Verzweigungpunkten, Math. Ann. bf 55 
(1902) 51-60; 

R. Vakil, Enumerative geometry of curves via degeneration methods, Harvard Ph.D. thesis (1997); 
I. Goulden and D. Jackson, Transitive factorisations into transpositions and holomorphic mappings on 
the sphere, Proc.Amer.Math.Soc. 125 (1997) 5160, math/9903094; 

S.Lando and D.Zvonkine, On multiplicitites of the Lyashko-Looijenga mapping on the discriminant 
strata, Funk. Anal. Appl. 33 3 (1999) 178-188; Counting ramified coverings and intersection theory on 
spaces of rational functions. I, math. AG/0303218; 

S.Natanzon and V.Turaev, A compactification of Hurwitz space, Topology, 38 (1999) 889-914; 
A.Okounkov, Toda equations for Hurwitz numbers, Math. Res. Lett. 7 (2000) 447-453; 
A.Givental, Gromov-Witten invariants and quantization of quadratic hamiltonians, math/0108100; 
S.Lando, Ramified coverings of the two-dimensional sphere and intersection theory in spaces of mero- 
morphic functions on algebraic curves, Russ.Math.Surv., 57 (2002) 463-533; 

A.Okounkov and R.Pandharipande, Gromov-Witten theory, Hurwitz theory, and completed cycles, 
Ann. of Math. 163 (2006) 517, math. AG/0204305; 

T.Graber and R. Vakil, Hodge integrals and Hirwitz numbers via virtual localization, Compositio Math., 
135 (2003) 25-36; 

M.Kazarian and S.Lando, Towards the intersection theory of Hurwitz spaces, math. AG/0410388; An 
algebro- geometric proof of Witten's conjecture, math/0601760; 

M.Kazarian, KP hierarchy for Hodge integrals, based on the talk at the Moscow Workshop on Combi- 
natorics of moduli spaces, Hurwitz numbers and cluster algebras (June 2008), http://www.mi.ras.ru/ 
~kazarian/papers / newwit0703 .pdf ; 

A.Mironov and A.Morozov, Virasoro constraints for Kontsevich-Hurwitz partition function, JHEP 
0902 (2009) 024, arXiv:0807.2843 

[15] J. Harer, D. Zagier, The Euler Characteristic of the Moduli Space of Curves I, Inv. Math. 85 (1986) 
457-485; 

S.K. Lando, A.K. Zvonkine, Graphs on Surfaces and Their Applications, Springer (2003); 
E.Akhmedov and Sh.Shakirov, Gluing of Surfaces with Polygonal Boundaries, to appear in Funkts. 
Anal. Prilozh., arXiv:0712.2448; 

A.Morozov and Sh.Shakirov, Harer-Zagier correlation functions in Gaussian matrix models , to appear 

[16] V.Dolotin and A.Morozov, Introduction to Non-Linear Algebra, World Scientific, 2007, hep- 
th/0609022; 

A. Morozov and M. Serbyn, Non-Linear Algebra and Bogolubov's Recursion, Thcor.Math.Phys. 154 

(2008) 270-293, hep-th/0703258; 



28 



V.Dolotin, A.Morozov and Sh.Shakirov, Aoc-structure on simplicial complexes, arXiv:0704.2609; 
Higher Nilpotent Analogues of Structure, Phys.Lctt. B651 (2007) 71-73, arXiv:0704.2884; 
A.Morozov and Sh.Shakirov, Analogue of the identity LogDet = TraceLog for resultants, 
arXiv:0804.4632; Resultants and Contour Integrals, arXiv:0807.4539; 

A.Anokhina, A.Morozov and Sh.Shakirov, Resultant as Determinant of Koszul Complex, 
arXiv:0812.5013 

[17] N.Amburg, E.Kreines and G.Shabat, Parasitic solutions of systems of equations determining the Belyi 
functions of plane trees, Mosc. Univ. Math.Bull. 59:1 (2004) 20-24; 

N.Adrianov, N.Amburg, V.Dremov, Yu.Levitskaya, E.Kreines, Yu.Kochetkov, V.Nasrctdinova, and 
G.Shabat, Catalog of dessins d'enfants with < 4 edges, arXiv:0710.2658 

[18] A.Morozov, String Theory, What is it?, Sov. Phys. Usp. 35 (1992) 671-714; 

A. Mironov, A. Morozov and L. Vinet, On a c-number quantum r-function, Teor.Mat.Fiz. 100 (1994) 
119-131 (Theor.Math.Phys. 100 (1995) 890-899), hcp-th/9312213; 

A.Gerasimov, S.Khoroshkin, D.Lebedev, A. Mironov and A.Morozov, Generalized Hirota Equations 
and Representation Theory. I. The case of SL(2) and SL q {2), Int.J.Mod.Phys. A10 (1995) 2589-2614, 
hep-th/9405011; 

S.Kharchev, A. Mironov and A.Morozov, Non-Standard KP Evolution and Quantum t -function, q- 
alg/9501013; 

A.Gorsky, I.Krichever, A.Marshakov, A. Mironov and A.Morozov, Integrability and Seiberg-Witten Ex- 
act Solution, Phys.Lctt. B355 (1995) 466-474, hcp-th/9505035 

[19] S.Kharchev, A.Marshakov, A. Mironov and A.Morozov, Landau- Ginzburg Topological Theories 
in the Framework of GKM and Equivalent Hierarchies, Mod.Phys.Lett. A8 (1993) 1047-1062, 
Theor.Math.Phys. 95 (1993) 571-582, hcp-th/9208046 

[20] A.Zamolodchikov, Theor.Math.Phys. 63 (1985) 1205; 

V.Fateev and A.Zamolodchikov, Nucl.Phys. B280 [FS18] (1987) 644; 

A.Gerasimov, A.Marshakov and A.Morozov, Hamiltonian Reduction of Wess-Zumino-Witten Theory 
from the Point of View of Bosonization, Phys. Lett. 236B (1990) 269-272; Free Field Representation 
of Parafermions and Related Coset Models, Nucl.Phys. B328 (1989) 664-676; 
A.Marshakov and A.Morozov, A Note on W 3 -Algebra, Nucl.Phys. B339 (1990) 79-94; 

A. Morozov, On the Concept of Universal W Algebra, Nucl.Phys. B357 (1991) 619-631 

[21] L. Chekhov and Yu.Makeenko, A Hint on the External Field Problem for Matrix Models, Phys. Lett. 
B278 (1992) 271-278, hep-th/9202006; The Multicntical Kontsevich-Penner Model, Mod.Phys.Lett. 
A7 (1992) 1223-1236, hcp-th/9201033; 

J. Ambjorn, C. Kristjansen and Yu. Makeenko, Higher Genus Correlators for the Complex Matrix 
Model, Mod.Phys.Lett. A7 (1992) 3187-3202, hep-th/9207020 

[22] B. Sagan, The Symmetric Group. Representations, Combinatorial Algorithms, and Symmetric Func- 
tions, 2nd edition, Springer- Verlag, New York (2001); 

B. Simon, Representations of Finite and Compact Groups, Graduate Studies in Mathematics, vol. 10, 
American Mathematical Society, Providence, RI (1996) 

[23] B.Eynard, All orders asymptotic expansion of large partitions, arXiv:0804.0381; 

A. Klemm and P. Sulkowski, Seiberg-Witten theory and matrix models, arXiv:0810.4944 

[24] V.Knizhnik, Multiloop amplitudes in the theory of quantum strings and complex geometry, 
Usp.Fiz.Nauk 159 (1989) 401-453 (Sov.Phys.Usp. 32 (1989) 945-971); 



29 



A.Gorsky, A.Marshakov, A.Mironov and A.Morozov, RG Equations from Whitham Hierarchy, 
Nucl.Phys. B527 (1998) 690-716, hcp-th/9802007; 

H.Bradcn, A.Mironov and A.Morozov, QCD, Wick's Theorem for KdV T-functions and the String 
Equation, Phys.Lett. B514 (2001) 293-298, hep-th/0105169 



30 



